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Abstract 


We calculate the contributions of penguin diagrams with internal u or c quarks to various inclu¬ 
sive charmless B-decay rates. Eurther we analyze the influence of the chromomagnetic dipole 
operator Qg on these rates. We find that the rates corresponding to 5 —> Xuus, B 
B —> Xsss, B —> Xssd and B X^^d dominated by the new penguin contributions. The 
contributions of Qg sizably diminish these rates. Despite of an increase of the total charmless 
decay rate by 36 % the new contributions are not large enough to explain the charm deficit ob¬ 
served by ARGUS and CEEO. We predict ric = 1.33 ± 0.06 for the average number of charmed 
particles per B-decay in the Standard Model. Then the hypothesis of an enhancement of the 
chromomagnetic dipole coefficient Cg by new physics contributions is analyzed. We perform a 
model independent fit of Cg to the experimental data. If the CKM structure of the new physics 
contribution is the same as in the Standard Model, \Cg{Myy) \ must be enhanced by a factor of 9 
to 16 in order to explain the observed charm deficit. 
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1 INTRODUCTION 


1. Introduction 

Precision measurements performed at the T (dS*) resonance find less charmed particles in the final 
states of B meson decays than theoretically expected. The CLEO 1.5, CLEO II and ARGUS data 
give [1] 

= 1.15 ±0.05 (1) 

for the average number of charm (anti-)quarks per B+/B°-decay. Complementary information 
on inclusive B decays can be obtained from the semileptonic branching ratio. The CEEO and 
ARGUS groups [1,2] have measured 

Bs^l = 10.23 ±0.39%. (2) 

The increasing experimental precision achieved in the current decade has been paralleled by a 
substantial progress in the theoretical understanding of the inclusive decay rates entering Bsl 
and Tie. Here the calculational key is the heavy quark expansion (HQE) [3, 4] of the decay 
rate in question: The leading term of the HQE reproduces the decay rate of a b-quark in the 
QCD-corrected parton model. The first non-perturbative corrections are suppressed by a factor 
of {I^QCD/rnhf and affect the rates by at most a few percent. Theoretically spectator effects of 
order lQT[‘^{kQCD/Tnbf [5,6] could be larger [6], but for the decay rates of B"^ and B° entering 
(2) and (1) they are experimentally known to be at the percent level as well [7]. The apparent 
smallness of these non-perturbative terms has shifted the focus towards the perturbative correc¬ 
tions to the free quark decay. The calculation of such short distance effects starts from an effective 
hamiltonian, whose generic form reads 

VCKM X! ^jQj ~ ^CKM fX! 

j=l \k=3 

Here is the Eermi constant and Vckm and products of elements of the Cabibbo- 

Kobayashi-Maskawa (CKM) matrix. The Wilson coefficients Cj encode the physics connected to 
the weak scale and play the role of effective coupling constants of the local interactions described 
by the operators Qj. Their precise form depends on the flavour structure of the decay and will be 
given below in (13). 

Decays with three different flavours in the final state such as b ^ cud can only proceed through 
the current-current operators Qi and Q 2 . T{b ^ cud) has been calculated to order as, which is 
the next-to-leading order (NEO), in [8]. The generic Eeynman diagram for these corrections are 
depicted in Eig. 1. In [9] the same diagrams have been calculated for r(5 — ccs) + T{b ^ ccd). 
The latter decays and the charmless nonleptonic decays, however, also involve penguin effects. 
Diagrams with insertions of the penguin operators have been taken into account only in 
the leading order (EO) [9-11], because their coefficients are much smaller than Gi 2 (cf. 
Tab. 1). The results for Bsl and ric read 

= (11.7 ± 1.4 ± 1.0)%, ne = 1.34^0.06. (4) 

Here the result for Bsl has been obtained in [10,12] with the analytical input from [8-10]. The 
second error bar has been added to account for the spectator effects estimated in [6]. Apparently 
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there is no speetaeular diserepaney between (2) and (4). The result for Uc in (4) does not involve 
the ealeulation of r(6 — ccd) + r(6 — ccs), but instead uses the experimental information on 
Bsl in (2) as proposed in [13,14]. We diseuss this in more detail in seet. 2. 

The diserepaney between (1) and (4) eonstitutes the “missing eharm puzzle”. The seareh for a 
theoretieal explanation has reeently foeused on new positive eontributions to the yet unmeasured 
eharmless deeay modes entering (4). Indeed, in a reeent analysis [14] V{h ^ no charm) has been 
estimated indireetly in two different ways: First the experimental information on final states with 
hadrons eontaining a c quark has been used and seeond data on deeay produets involving a c 
quark have been analyzed. For the CLEO data both methods eonsistently indieate an enhanee- 
ment of r(6 —no charm) by roughly a faetor of 14 eompared to the theoretieal predietion 
in [11]. 

Next we briefly diseuss the LEP results for Bsl and ric [15]. The EEP Z-peak experiments 
eneounter a mixture of b-flavoured hadrons. In order to allow for a eomparison with (2) one must 
eorrect the EEP result Bsl^^ = 10.95 ± 0.42% [1] for the different lifetimes [7] of the hadrons 
in the mixture [16]: 

^z,corr,exp ^ 11.13±0.42%, = 1.22 ± 0.08. (5) 

These data are eonsistent with the theory (of. (9) below), but the analysis in [14] has found 
evidenee for an enhanoed V{h ^ no charm) also from the EEP data. Eurther the two methods 
used in [14] have given less eonsistent results for the EEP data than for the CEEO data. In 
addition the EEP measurements involve the baryon, whose lifetime is either not properly 
understood theoretieally or inoorreetly measured. (If the latter is the ease, in (5) must 

be replaeed by the unoorreoted Bgi^^, whieh reduees the 2(j diserepaney between (2) and (5).) 
Henee in our analysis we will mainly use (2) and (1). 

Now two possible souroes of an enhanoed V{h ^ no charm) are ourrently disoussed: The authors 
of [17] stress the possibility that the Wilson coeffioient Cs of the ohromomagnetie dipole operator 
Qs is enhanoed by new physies eontributions. On the other hand in [14,18] an explanation within 
QCD dynamics is suggested: An originally produced (c, c)-pair can annihilate and thereby lead 
to a charmless final state. 

The calculation of matrix elements involving Qs-e and Qs does not exhaust all possible penguin 
effects. In this paper we calculate the contributions of penguin diagrams with insertions of the 
current-current operator Q 2 to the decay rates into charmless final states (see Pig. 2). Such a 
penguin diagram with a (c, c)-pair in the intermediate state involves the large coefficient C 2 
and the CKM factor Vcb 3> \Vub\- It is precisely the short distance analogue of the mechanism 
proposed in [14,18] and surprisingly has not been considered in the perturbative calculations 
[9-11,20] of the decay rates entering Bsl and Uc- 

Eurther we calculate the diagrams involving the interference of the tree diagram with Qs in 
Pig. 3 with Qi_ 6, which belongs to the order as as well. The consideration of these diagrams is 
mandatory, if one wants to estimate the effect of an enhanced coefficient Cg on r(6 —no charm) 
proposed in [17]. 

The paper is organized as follows: In the following section we set up our notations and collect 
results from earlier work. The calculation of the penguin diagram contributions and the Qi^^-Qg 
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2 PRELIMINARIES 




Eigure I: Example of a 
current-current diagram. The 
cross denotes the inserted 
operator, d' equals d or s. 


Eigure 2: Penguin diagram 
involving Q 2 . The internal 
quark q' can be u or c. The 
corresponding diagram with 
Qi vanishes. 


Eigure 3: Tree diagram con¬ 
tribution of Qs to r(& ^ 
qqd'). 


interference terms is presented in sect. 3. The phenomenologically interested reader is refered to 
sect. 4, in which we discuss our numerical results. In sect. 4 we also comment on the mechanism 
proposed in [18]. Eurther a potential enhancement of Cg is analyzed by a model independent fit 
of Cg to the experimental data. Einally we conclude. 

2. Preliminaries 

2.1. Bsl and iic 

Eor the theoretical description of the various decay rates it is advantageous to normalize them to 
the well-understood semileptonic decay rate [13,14]: 



This eliminates the factor of m\ common to all decay rates. Eor the charmless decays we will 
further use 


r4 


E 

q=d,s 

q'=u,d,s 


^q'q'q 



T 2r^e 


(7) 


The semileptonic branching ratio reads 


Bsl 


1 

^ CT + + E, =d,s \T cuq ^ccg] 


( 8 ) 


Small contributions such as Tut = 0.004 and radiative decay modes have been omitted in (7). In 
(8) also Lucs = 0.05 [11] has been neglected. The numerical value for Bsl in (4) is the average 
of the two results given in [10] corresponding to two different renormalization schemes (on-shell 
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vs. MS [21] quark masses). In the corresponding expression for ric the term Vccs + rccd? which 
suffers from sizeable theoretical uncertainties, is eliminated in favour of Bsl [13,14]: 

2 ^2 -f -f "1“ f'cus “f 2^^^ Bgj^. (9) 

Yet with [8,22] 


Tens + Tend = 4.0 ± 0.4, Ter = 0.25 (10) 

and B'^gl in (2) one obtains 

= 1.36^0.04 - 0.205-r^. (11) 

Now the current-current type radiative corrections to (cf. Fig. 1) have been calculated in [8, 
20,23] using different renormalization schemes. The penguin operators Qz-g have been included 
within the LO in [11]. With up-to-date values for Bsl and the CKM elements the calculation 
of [11] yields = 0.11 ± 0.08. Inserting this result into (11) yields the numerical value in (4). 
In contrast the indirect experimental determination in [14] has found = 1.6 ± 0.4. In order to 
reproduce the experimental value of rze in (1) one needs = 1.0 ± 0.4. 

So far the penguin diagrams of Fig. 2 have not been calculated for all possible charmless B decay 
modes. Yet for the pure penguin induced decays of a 5-quark into three down-type (anti)-quarks 
the contribution of the diagram of Fig. 2 has been obtained in terms of a twofold integral repre¬ 
sentation in [19]. Likewise the effect of penguin diagrams on the analysis of CP asymmetries has 
been studied in [19,24] and in [25] penguin effects on exclusive decays have been studied. 


2.2. Decay rates to order 


In order to describe decays of the type b qqd one needs the following hamiltonian: 


H 


Gf \ 




C E 

j&v 




( 12 ) 


Here ■Cn + Cc + 6 = 0 due to the unitarity of the CKM matrix and V = {3,.. .6, 8}. H in (12) 
comprises the following operator basis 


Qi = {dq)v-A ■ {qb)v-A ■ 1 
Ql = {dq)v-A ■ {qb)v-A ■ 1 
Ql = {db)v-A ■ {qq)v-A ■ 1 , 
Ql = idb)v-A ■ iqq)v+A ■ 1 , 

Qj = E Q' 


with q = uor q = c, 

Ql = {db)v-A ■ {qq)v-A ■ 5 
Ql = idb)v-A ■ iqq)v+A ■ 2 
for 3 < j < 6, 


q=u,d,s,c,b 

Qs = da^" [rridL + mi,R] T% ■ G% . (13) 

^The overall sign of the matrix element {Qs) depends on the chosen sign of the coupling g in the covariant 
derivative in the QCD lagrangian. The definition in (13) complies with the result in [26,27], if the covariant deriva¬ 
tive is chosen as = d^ — igT°‘A^, so that the Feynman rule for the fermion-gluon vertex is igT°‘. By convention 
the notation Qi is reserved for the magnetic (i& 7 -operator, which we do not need in our calculation. 
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2 PRELIMINARIES 


The colour singlet and non-singlet structure are indicated by 1 and 1 and V ± A is the Dirac 
structure, i.e. 


(1-75) (1-75) &«• 

Next it is useful to expand the renormalized matrix elements in Ug and to separate the result from 
current-current diagrams (see Eig. 1) and penguin diagrams (see Eig. 2): 


(OTdigfit) = {g|')<"' + j^[(Q(>w + {g»')li>J+o(aj), 3 = 1 , 2 , (w) 


, ^ j I peng 


^jk (Qfc)®, 


p = Pb-Pd 


(15) 


ker 


Of course ( Q'j and ( are non-zero only for q = q' = u, recall that we do not consider 

g = c in this work. In (15) we have expressed the result of the penguin diagram in terms of 
the tree-level matrix elements. There n is the renormalization scale. Eor the momentum flow 
cf. Eig. 2. 

The quark decay rate is related to the matrix element of H via 

T{b^qqd) = ^ j ^ ^ (5^^^ {pb + Pq -Pq- Pd) ( qqd\H\b){ qqd\H\ b)*. 

2mbJ {2TTf8\EqEqEd\ 

The bar over ( qqd \ H\b){ qqd \ H\b)* denotes the average over initial state polarizations and the 
sum over final state polarizations. Next we expand the decay rate to order 

T{b^qqd) = + ^^[ATg, + AT,g^g + ATs + ATw]+0{al). (16) 

Eor b qqs decays one simply substitutes d by s in (13-16). Now the first two terms of the 
NLO correction in (16) describe the effect of current-current and penguin diagrams involving 
Qi or Q 2 - ^rs likewise contains the matrix elements of Qs- The remaining part ATpy of the 
NLO contribution is made of the corrections to the Wilson coefficients [20,28] multiplying the 
tree-level amplitudes in We write 

Ciifi) = + J = l,...6. (17) 

Here ACj is the NLO correction to the Wilson coefficient. ACj depends on the renormalization 
scheme chosen. This scheme dependence cancels with a corresponding one in the results of the 
loop diagrams contained in ATcc and ATp^ng- Eor example the scheme dependence of ACy 2 
cancels in combination with the current-current type corrections to Qi and Q 2 of Eig. 1. Since 
we do not include the unknown radiative corrections to the penguin operators Qs-e in (16), we 
must likewise leave out terms in ACj related to the NLO penguin-penguin mixing in order to 
render T in (16) scheme independent. We ban these technical details into the appendix. The 
values of the Wilson coefficients needed for the numerical evaluation of the various decay rates 
are listed in Tab. 1. 
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j 

1 

2 

3 

4 

5 

6 

8 

(/i = mb) 

= mb) 

m 

1.1077 

1.0731 

|an||| 


Blil 

nn»y 

-0.1495 

Cf\ix = mb/ 2) 

= mb/2) 


m 

M 


0.0100 

0.0149 


-0.1663 

Cf\n, = 2 mb) 

= 2 mb) 

m 

IH 

0.0071 

0.0073 

nxo^ 

0.0054 

0.0079 

nxig«j 

-0.1355 

ACjifi = mb) 

2.719 

-1.744 

0.380 

-0.1050 

-0.223 

0.384 



Table 1: Wilson coefficients used in our analysis. Cf'* is the LO expression, is the NLO 

coefficient in the NDR scheme. In above the NLO corrections to penguin-penguin mixing 

have been omitted in order to render V in (28) scheme independent as described in the text. For 
jj, = mb this affects C 3 and C 5 by 12 % and 25 %, but is negligible for the other coefficients. 
cf'* and are needed for the numerical evaluation of the decay rate in (28). ACj in the last 

line is defined in (42). The top and bottom mass are chosen as mt = = 168 GeV and 

mb = 4.8 GeV. Further as{Mz) = 0.118 [29], which corresponds to 0 : 5 ( 4 . 8 GeV) = 0.216. In 
the table is the scheme independent coefficient mentioned in the appendix. 


In the LO the decays b sss, b ssd, b dds and b ddd can only proceed via Qs-e and 
Qs, while b uud and b uus also receive contributions from Qi and Q 2 . We combine both 
cases in 


p(o) ^ Glml 

6471-3 


2 


6 


12 . \ucrcf%, + ^ 




i,j=l 


ij=3 


- 2 t ^ Cf’cf Re (C6) 


i=l,2 
7=3,...6 


with t = 1 for g = M and t = 0 for g = d, s. The b^’s read 


bij 


with Qi 2 = Qi 2 here- Setting the final state quark masses to zero one finds 


bij = 


1 + r/3 for i, j < 4, and i + j even 
1/3 + r for i, J < 4, and i + j odd 


^55 — ^66 
^56 = ^65 


(18) 




= 1 , 

= 1/3. 


( 20 ) 


Here r = 1 for the decays b ddd and b —> sss, in which the final state contains two identical 
particles, and r = 0 otherwise. The remaining bi/s are zero. Clearly for the g 7 ^ m the bi/s as 
defined in (19) vanish, if f < 2 or j < 2. Yet in the formulae for the decay rate we prefer to stress 

































8 


2 PRELIMINARIES 



Eigure 4: The two diagrams contributing to in (18). The crosses represent any of Qi-e- 


this fact by keeping the parameter t, which switches the current-current effects off in the penguin 
induced decays. Our results in (20) agree with the zero mass limit of [9,10]. The bi/s of (19) are 
visualized in Eig. 4. 

AT w simply reads 


AT^^ — 


G^ml 

647r^ 


t E 


*,i=i 


C^ACj 


t E Re (C6) h, 


+ E 


i,j=3 


( 0 ) 


C^ACi 


bij 


i=l,2 

i=3,...6 


( 21 ) 


The current-current type corrections proportional to C [^2 ' [8,20,23] 


AT. 


= t 


Gpml 

6471-3 


2 le. 


2 

’E 

*4=1 


’^ij 


( 22 ) 


with t defined after (18). In the NDR scheme with the standard definition of the evanescent 
operators [30] the diagrams of Eig. 1, the bremsstrahlung diagrams and the subsequent phase 
space integrations yield [8]: 


j^NDR 

'111 




'22 





'12 





'21 




17 

y 



The inclusion of ATcc is necessary to fix the definiton of the b-quark mass entering in (18). 
The values quoted in (23) correspond to the use of the (one-loop) pole quark mass in (18). 

In the same way we write 


AT 


peng 


Gpml 

647r^ 


2 Re 


t Y, Cfcf {.|&.;fe)+C9.i(0)l 


- cfcf&[e9„(ij+e:9«(o)] 

i=l,2 
7=3,...6 


( 24 ) 
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/ / 

Figure 5: The diagrams contributing to ATpeng in (24). The left cross denotes Ql or with 
q' = u,c and the right cross represents any of Qi-e- The dashed line indicates the final state with 
q being u, d or s. 



Here gij is visualized in Fig. 5 and defined by 


9ij 





d^Pqd^Pqd^Pd 

{2TT)^8\EqEqEd 


( 5 ^^) {Pb + Pq 


Pq - Pd) ■ { Qi * (25) 


with Xq/ = TUq'/mb- We do not display the /i-dependence of the Cj's, hi/s and Pi/s in formulae 
for the decay rate such as (22) or (24) to simplify the notation. Now ATpeng in (24) is more com¬ 
plicated than AFcc in (22) for two reasons: First interference terms of different CKM structures 
appear and second the internal quark in the penguin graph of Fig. 2 can be g' = c or g' = u. 
Further the charm quark mass enters Pij with Xc = rric/mb. 

Finally AFg in (16) is given by 


Ar« = 


G%m 


647r^ 


2 Re 


-*a.crEcf'>.s + ie 

i=i 


2 ^(0) 

t\ ^8 9j8 

i=3 


(26) 


Here the tree-level diagrams with Qs already contribute to order 

The phase space integrations are contained in the coefficients bjg in (26). They are depicted in 
Fig. 6 and are defined as 




IGtt® 47r 
mf as 


d^Pqd^Pqd^Pd 

{27r)^8\EqEqEd 


5 ^^) {pb +Pq-Pq- Pd) * 


bsj. (27) 


It is instructive to insert the above expressions for AFcc, ATpeng and AFg into (16). The decay 
rate then reads 


r = 


Gpuil 

647r® 


Re h ^ G,Gq 


*j=i 


l 2 r , (p) 1 ^ |2 


I'Cm I hj + 


47r 

(/^) 


|C«I 2 [hij -f gij{0) 9ij{^c)] 

2 9ij i^c 


47r 
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Eigure 6: The diagrams contributing to AEg in (26). The left cross denotes Qg and the right cross 
represents any of Qi-e- 


2 ctc, 


i=l,2 

i=3,...6 


«a. K +ife(o) - mM] - i&P 


An 


An 


Y. Cfij 161 " h, 

i,j=3 

(/^) 


An 


-Cg 


2 6 
|2 


«e6i:c,2(,,s+i6ri:cj2(,,s 

j=i i=3 


(28) 


Here the Cfs are the Wilson coefficients of (17) including NLO corrections. Cf^ rather than 
CNDr sj^ould be used in the terms of order as in (28) for consistency. The unitarity of the CKM 
matrix has been used to eliminate from (28). 

Erom (24) or (28) one notices that penguin induced decays with f = 0 receive radiative correc¬ 
tions proportional to the large coefficient C 2 , which does not enter the tree-level decay rate in 
(18). Eurther ATp^ng depends on the CKM phase 6 , because 6 and the loop functions Qij are 
complex. Decay rates and CP asymmetries for these penguin induced decays have been derived 
in [19] in terms of a two-fold integral representation taking into account the interference of the 
penguin diagram involving Q 2 with Qi-e- In sect. 3 we derive analytical results for the decay 
rates and also include ATg. 

In the decays b uud and b —> uus the main focus is on the first sum in ATpeng in (24) con¬ 
taining products of two of the large coefficients Ci and C 2 . We also keep the second sum in (24), 
but remark here that these terms are not the full set of one-loop radiative corrections involving 
one large coefficient Ci ^2 and one small coefficient C 3 -Q: In addition to the radiative corrections 
calculated in this paper there are also current-current diagrams (see Pig. 1) and penguin diagrams 
(see Pig. 2) with penguin operators Qg-e- In decays with a (m, M)-pair in the final state these ma¬ 
trix elements interfere with those of and Q 2 and therefore also yield a term proportional to 
<3*1,2 ■ C's-e- 

The terms proportional to Cj ■ Cg comprised in AFg are interesting in order to confirm or falsify 
the mechanism proposed in [17]. If new physics indeed dominates Cg, then T{b s + g) oc \Cg\^ 
considered in [17] is enhanced. Yet AFg is linear in Cg, so that the sign of the new physics 
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contribution may determine whether r(b ^ no charm) is enhanced or diminished. A calculation 
similar to ours has been partly done in [31]. Yet in [31] some questionable approximations have 
been made: For example the operator mixing has been neglected and the top quark has not been 
integrated out but instead formally treated as a light particle. In some cases the results of [31] 
differ substantially from those in [19]. 

We close this section with the formula relating V in (28) to the r^^^gg’s defined in ( 6 ). To this end 
we need the semileptonic decay rate to order as [32]: 

1 hsL^Xc) + O . 

The tree-level phase space function is 

/i(a) = 1 — 8a — 12a^ Ina + 8a^ — a^. 

The analytic expression for hsiixc) can be found in [8,32]. The approximation 
hsL (Xc) = -3.341 + 4.05 {xc - 0.3) - 4.3 {xc - 0.3)^ 


T{b 




ml 


cevp = 


1927r3 


\Vcb? h (x) 


holds to an accuracy of 1 permille in the range 0.2 < Xc < 0.4. Here Xc = nic/mb is the ratio of 
the one-loop pole masses. We further include the hadronic corrections to the free quark decay of 
order 1/m^ obtained from the HQE [3]. This yields 


f'qqd 


192%^ 


Glml iKfel fi {xD 


r + r® 


. hiXc) 


- 1 


mt 


asih) 

271 


hsL {Xc) 


+ (5r|.(29) 


Here A 2 = 0.12 GeV^ encodes the chromomagnetic interaction of the b-quark with the light de¬ 
grees of freedom. Further corrections are contained in (iF. It is obtained from by substituting 
bij with 6bij in the definition (18). One finds from [3]: 


6bi2 — Sbi4 — 5621 — bb23 — 6bs2 — <^^34 — bb^i — 5643 



-0.042 


for mb = 4.8 GeV, 


while Sb^Q = bb^^ is unknown yet. If there are no identical quarks in the final state, the remaining 
5bi/s vanish. Otherwise 

bbs 3 = 5644 = 5634 , ( 5&55 = bbee 7 ^ 0 for 6 —ddd and b sss (30) 


with bb^Q unknown. The dependence on Ai parameterizing the effect of the b-quark’s Fermi 
motion on the decay rates cancels in the ratio in (29). 

For the calculation of in (7) we also need Vue and Vgg. The corresponding expressions are 
[3,27]: 


r 


ue 


r 


sg 


Vcb\ 

VtlVts 

Vcb 


1 

77^!) r 

^ Sasjfi) 

fi{xl) 


[/iSL(O) - hsL{Xc)] + 6 
Zti 

[Csih)]"- 


fi{xl) 


A 2 

ml 
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3 CALCULATION 


3. Calculation 

This section is devoted to the calculation of the gi/s and fe^g’s entering (24,26,28). All results 
correspond to the NDR scheme. 

The first step is the same for all b-decays under consideration: The penguin diagram of Fig. 2 
must be calculated to obtain the Vi/s in (15): 



r26 = r23 = rss = -2^24- (32) 

Here p = pt — Pd the momentum flowing through the gluon leg and m is the internal quark 
mass. The infinitesimal “i(5”-prescription yields the correct sign of the imaginary part of the 
logarithm in the case p^ > Am? and likewise regulates the square root for p^ < Am^. The ri/s 
and r 28 are zero. 

Next we combine (25) and (15) to obtain the coefficients gij in (24) and (28): 


Qij I 


d^Pqd^Pqd^Pd .U) t , N 

{2.mEqE-qEdr ^^^^-^‘^-^^^ 


({Pb-Pdf,Xcmb,p) ( 33 ) 


The corresponding expression for an internal u-quark is obtained by substituting c with u in (33). 
For the decays b uud, b uus, b —> ssd and b dds one finds: 
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922 0 , 


/i 


921 X, 


mb, 
j^\ 


— 924, ( 0 , 


mtJ 


iL 

mb 


— 926 ( 0 , 


J^\ 

mb) 


mb 




mb 


3 4 mfe 2. 

-— In- m 

2 3 /i 3 




= 9211 [x, — = 92b [x, — = giAx, — =0, j = 1,... 6, (34) 


mb 


where Xc = mc/mb- Numerically one finds for actual quark masses: 

^ 22 ( 0 , 1 ) = -0.67- 0.93 f, ^ 22 ( 0 . 3 , 1 ) = -0.69 - 0.23 f. 


(35) 


The near equality of the real parts in (35) is a numerical accident. 

In the case of two identical particles in the final state both diagrams of Fig. 5 contribute. Then 
g 23 is no more zero, but 


/i 


mb 


923 \X, - = 922 X, — 


/r 


mb 


for b ddd and b sss. 


(36) 


The remaining gij's are as in (34). As an analytical check we have confirmed that the /i-depen- 
dence in (34) and (36) cancels with the /i-dependence in the Wilson coefficients in (18) to order 

We now turn to the calculation of ATg in (26). Performing the phase space integration in (27) 
yields 

b - b - b - 

^18 = ^38 = ^58 = 0 (37) 

ssd and b dds. If identical particles are present, 633 is no more 


for b uud, b uus, b 
zero, but instead reads 


hs — 


16 

y 


for b ddd and b sss. 


(38) 


The remaining bjgs are as in (37). 


4. Charmless decay rates 


4.1. Standard Model 


In this section we discuss our numerical results for the various decay rates. We use the following 
set of input parameters: 


Vcb 


0.08 ± 0 . 02 , 


a, (Mz) = 0.118, 


5 = 90° ± 30°, Xc = 0.29 ± 0.04, ^ = mb = 4.8 GeV, 

\Vcb\ = 0.038, mt (mt) = 168 GeV. (39) 


The rqqqds sizeably depend on \ Vub/Vcb\, d, Xc and especially on the renormalization scale, which 
will be varied in the range mb/2 < gi < 2mb. The quark masses in (39) are taken from [33], the 
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4 CHARMLESS DECAY RATES 



input 

uud 

uus 

dds 

fir 

sss 

lal state 

ssd 

ddd 

no eharm 

'Tqqq' 

as in (39) 

0.040 

0.021 

0.018 

0.015 

8.9 ■ 10-^ 

7.2 ■ 10“^ 

0.14 

jj, = mb/2 

0.044 

0.033 

0.029 

0.024 

14.0 ■ 10“^ 

11.4-10-^ 

0.19 

fi = 2 mb 

0.037 

0.014 

0.011 

0.009 

5.5 ■ 10-^ 

o 

1 

4^ 

0.11 


Vub 

= 0.06 

0.023 

0.020 

0.018 

0.015 

OO 

o 

1 

4^ 

o 

1 

4^ 

0.11 


Lffr 

Lb 

= 0.10 

0.062 

0.023 

0.018 

0.015 

o 

1 

4^ 

7.5 ■ 10-^ 

0.18 

5 = 60° 

0.044 

0.017 

0.018 

0.015 

5.9 ■ 10"^ 

4.8 ■ 10“^ 

0.14 

5 = 120° 

0.036 

0.025 

0.018 

0.014 

12.2 ■ 10“^ 

o 

o 

o 

1 

4^ 

0.14 

Xe = 0.25 

0.034 

0.019 

0.016 

0.013 

OO 

o 

1 

o 

1 

4^ 

0.12 

Xc = 0.33 

0.048 

0.023 

0.020 

0.017 

o 

1 

4i- 

o 

1 

4^ 

0.16 

Br 

as in (39) 

0.41% 

0.22% 

0.18% 

0.15% 

o 

1 

CO 

7.4 ■ 10-3% 

1.4% 


Table 2: The values of rqqq/ for the various final states as defined in (6). The input parameters are 
ehosen as in (39) exeept for the quantity listed in the seeond eolumn. The last eolumn lists 
defined in (7). Br in the last row is the branehing ratio for B ^qqq' obtained by multiplying 
Tqqq/ With BsL = 10.23%. 


values for \Vub/Vcb\ and \ Vcb\ have been presented in [34]. The range for the CKM phase 6 has 
been obtained from the NEO analysis of €k and Auib in [35]. The dependenee on as {Mz) in 
the range 0.112 < o;* {Mz) < 0.124 [29] is weaker than the /i-dependenee. Our results are listed 
in Tab. 2. 

Keeping the physieal input parameters as in (39) and varying the seale in the range mf,/2 < /i < 
2mb the eharmless deeay modes sum to 


= 0.15 ±0.04. 

The values for Vgg and Vue entering this result are 

Tsg = 0.02 ± 0 . 01 , r^e = 0.01 ± 0 . 00 . 

Ineorporating also the uneertainties in \ Vub/Vcb\ and Xc one finds 


= 


0.15 ±0.08. 


(40) 
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final state 

^(0) ^(0) 
CX 0^ 2 ^1,2 

^(0) ^(0) 

OC <-/3_6 V/1_0 

(X. AT 

oc ATcc 

CX pQjig 

oc ATg 

uud 

107 

-4 

-5 

6 

-5 

1 

VMS 

10 

39 

-7 

1 

72 

-14 

dds 

0 

46 

-8 

0 

78 

-16 

sss 

0 

44 

-20 

0 

92 

-16 

ssd 

0 

54 

-9 

0 

74 

-19 

ddd 

0 

52 

-24 

0 

91 

-19 

sg 

0 

0 

0 

0 

0 

100 

no charm 

46 

17 

-6 

2 

31 

11 


Table 3: Separate contributions to r(6 —> percent of the total rate for the input parameters 

of (39). In the third column the contribution from the LO matrix elements of penguin operators 
and their interference with matrix elements of (^ 1,2 is shown. In the last row r(6 — ueug) 
entering has been assigned to the second column listing the current-current part. 


To discuss the results for the individual charmless decay modes it is instructive to look at the 
separate scheme-independent contributions from T^^^ ATcc, ATpeng, ATg and AT w (cf. the ap¬ 
pendix) to the decay rate. These contributions are listed in Tab. 3, in which also the contributions 
from penguin operators to are shown. 

All decays except for B X^ud are dominated by penguin effects. The sizeable /i-dependence 
in these decays can be reduced by calculating the current-current type radiative corrections to 
penguin operators. Note that all these decay rates are even dominated by the penguin diagrams 
calculated in this work. On the other hand the terms stemming from Qg lower the penguin in¬ 
duced rates. InB ^ no charm, however, the net effect of Qs is positive because of the two-body 
decay b ^ sg proportional to C|. 

The calculation of in [11] has included T^'^^ AT w, and ATcc- Taking into account that in [11] 
a high (theoretical) value for Bsl has been used, the result of [11] translates into 


r4 


0.14 


for 




0 . 10 . 


The corresponding value in Tab. 2 is = 0.18 showing the increase due to ATpeng + ATg. 
Despite of the 36 % increase in in (40) the value is still much below = 1.0 ± 0.4 needed to 
solve the missing charm puzzle. The new theoretical prediction for ric is 


n, 


1.33^0.06, 
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4 CHARMLESS DECAY RATES 


which is only marginally lower than the old result in (4). Still an enhaneement of by a faetor 
between 2.6 and 20 is required. Yet with the result in (34) for the penguin diagram of Eig. 5 at 
hand we ean estimate the non-perturbative enhaneement of the (c, c) intermediate state neees- 
sary to reproduee the experimental result. The meehanism proposed in [14,18] eorresponds to a 
violation of quark-hadron duality, whieh we may parametrize by multiplying the (c, c) -penguin 
funetion gij {xc, /i/rrib) in (28) by an arbitrary faetor d. We find that d must be ehosen as large as 
20 in order to reaeh the lowest desired value = 0.6 for the eentral set of the input parameters 
in (39). Yet in this ease one must also inelude the double-penguin diagram obtained by squar¬ 
ing the result of Eig. 2. Taking this into account non-perturbative effeets must still inerease the 
(c, c) -penguin diagram by roughly a factor of 9 over its short distance result in the NDR seheme. 
Having in mind that the phase spaee integration eontained in {xc, fi/rri},) implies a smearing 
of the invariant hadronie mass of the (c, c)-pair sueh a large deviation of quark-hadron duality 
seems unlikely. 

Our results for the branehing fraetions of & — dds and b sss differ by roughly a factor of 
1/2 from the results in [19]. The main souree of the discrepaney is the different ealeulation of 
the total rate Ttot entering the theoretical definition of the branching fraction. Our predictions 
for the branehing ratios in Tab. 2 contain the total rate via Ttot = Tsl/B^sl^ while in [19] 
^tot = EcftP/(647r^) has been used. This approximation negleeting the RG effeets appears 

to be too erude and is responsible for 50 % of the diserepaney. The remaining differenee is due to 
the effeet of ATg not eonsidered in [19] and the use of different values of the Wilson eoeffieients. 
Our predietions for Br{h ssd) and Br{h ddd) are even smaller by a faetor of 5 than the 
results of [19]. This is due to the faet that in addition the value of \ Vtd/Vcb\‘^ used in [19] is more 
than twiee as large as the present day value used in our analysis. We remark that with our new 
results the number of B-mesons to be produced in order to deteet the inelusive CP asymmetries 
eorresponding to these deeays is substantially larger than estimated in [19]. 


4.2. New physics 

In the Standard Model the initial eonditions for Cg-e and Cg are generated at a seale fx = 0{Mw) 
by the one-loop bsg-\ertex funetion with a kE-boson and a top quark as internal partieles. Due 
to the helieity strueture of the eouplings positive powers of rrit are absent in the bsg-\ertex. In 
many extensions of the Standard model sueh a helieity suppression does not oeeur [17]. Eurther 
in supersymmetrie extensions flavour ehanging transitions ean be mediated by gluinos, whose 
eoupling to (s)quarks involves the strong rather the weak eoupling eonstant. Reeently a possi¬ 
ble enhancement of Cg affeeting r-gg and thereby has attraeted mueh attention [17]. In the 
following we will perform a model independent analysis of this hypothesis. 

New physies affeets the initial eondition of the Wilson eoeffieients ealeulated at the seale of the 
masses mediating the flavour-ehanging transition of interest. We will take fx = Mw as the initial 
seale for both the standard and non-standard eontributions to the Wilson eoeffieients. The renor¬ 
malization group evolution down to jx = rrib mixes the various initial values and ean damp or en¬ 
hance the new physies effeets in Ci{Mw )-Eor example the Standard Model value of Cg(4.8 GeV) 
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ACgCM^) 


Figure 7: vs. ACs{Mw) parameterizing new physics contributions to C^{Mw)- The dark 

shading marks the region = 1.0 ± 0.4 needed to reproduce the experimental result for Uc in 
(1). The lightly shaded area corresponds to = 1.6 ± 0.4 obtained from the analysis in [14]. 


in Tab. 1 is mainly a linear combination of C 2 {M]y) and C^{Mw): 

Csirrib) = —0.08(7*2 

with Cs{Mw) = —0.1. If one enlarges Cs{Mw) by a factor of 10 while keeping Ci_q{Miv) 
fixed, grows by a factor of 5.5. The sensitivity of to C 3 -e{Mw) is much smaller, increasing 
the latter by a factor of 10 enhances only by a factor of 1.6. Hence in the following we will only 
focus on (7*8 {Mw) = —0.1 + ACs{Mw), where ACs{Mw) is the new physics contribution. For 
simplicity we will further assume that the CKM structure of the new contributions is the same 
as in the Standard model and neglect the possibility of new CP-violating phases by assuming 
AC^{Mw) to be real. 

In Fig. 7 is plotted versus AC^{Mw). Solving for = 1.0 ± 0.4 yields two solutions: 


AC,{Mw) = - 1.2 , 


ACsiMw) = 1.5 ±0.3. 


(41) 
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A SCHEME INDEPENDENCE 


The central values correspond to an enhancement of the Standard Model value for Cs{Mw) 
by factors of 13 and (—14). We hope to resolve the twofold ambiguity after calculating the 
contribution of Qs toT{B ^ Xccs)- 
The EEP data in (5) correspond to = 0.4 ± 0.5 and 

ACsiMw) = -0.5 , ACsiM^) = 0.9 +;;;6 

showing the consistency of (5) with the Standard Model. 

5. Conclusions 

We have calculated two new contributions to the inclusive decay rates of B-mesons into vari¬ 
ous charmless final states: Eirst we have obtained the results of penguin diagrams involving the 
operator Q 2 and a c- or w-quark in the loop putting special care on the renormalization scheme 
independence of our results. Second we have calculated the influence of the chromomagnetic 
dipole operator Qs on these decays. The former contributions have been found to dominate the 
branching fractions for B ^ Xuus, B ^ X^js^ B Xsss, B ^ Xssd and B X^2d- 
effect of Qs on these decay modes is also sizeable and decreases the rates. On the other hand the 
decay rate for B X^ud is only affected by a few percent. Our results increase the theoretical 
prediction for Br{B no charm) by 36 %, which is not sufficient to explain the charm deficit 
observed in B-decays by ARGUS and CEEO. If a breakdown of quark-hadron duality due to in¬ 
termediate (c, c) resonances is to explain the “missing charm puzzle”, the phase space integrated 
penguin diagram with an internal c-quark must be larger than the perturbative result in the NDR 
scheme by roughly a factor of 9. 

We have then analyzed the hypothesis that new physics effects enhance the coefficient Cs of Qs 
and have performed a model independent fit of Cs to the experimental data on ric- The renor¬ 
malization group evolution from /i = Mw to jj, = rrib has been properly taken into account. 
One finds two solutions for Cs'. Eor the central values of the theoretical input and the data of 
the T(4S') experiments Cs{Mw) must be larger by a factor of 13 or (—14) than in the Standard 
Model, if the new physics contributions have the same CKM structure as the Standard Model 
penguin diagram. 
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A. Scheme independence 

Now we discuss the cancellation of scheme dependent terms between the NEO Wilson coef¬ 
ficients and the loop diagrams contained in ATcc and ATpeng- We define scheme independent 
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combinations AFcc, AFpeng and AF^y, which allow a meaningful discussion of the numerical 
sizes of these separate contributions to F as performed in sect. 4. Finally we comment on the 
scheme independence of AFg. 

The NLO correction to the Wilson coefficients in (17) can be split into two parts: 

AQ(/^) = + 7 = 1,... 6. (42) 

k=l 

ACj{iJ,) contains the contributions stemming from the weak scale. It is independent of the renor¬ 
malization scheme and proportional to as{Mw)/as{n). Yet the Jj^s in (42) are scheme depen¬ 
dent. The precise definitions of the terms in (42) can be found in [28, 30]. We now absorb the 
terms involving Jjk into AFcc and AFpeng, so that the latter become scheme independent. 

The identification of scheme independent combinations of one-loop matrix elements and Jj^’s is 
most easily done, if one expresses the loop diagrams in terms of the tree-level matrix elements. 
The combination 


(/, "ig', ft) + Jkj, 7 < 2 and /c > 3, (43) 

a' 

of the coefficients in (15) and the J^j’s is scheme independent [28]. Substituting with (43) in 
(33) one finds the scheme independent quantity: 


AF 


peng 


= AF 


peng 


+ 


Glml 

6471-3 


2 Re 


-i«.6 i: cf>cfjub,k 


fc=3,...6 


E cf’cfJid'ii, 


i=l,2 

j,k=3,...6 


(44) 


used. In the NDR scheme the J^/s in (44) evaluate to 

31 = 

-0.877, 

jNDR 

'-'32 ~ 

-0.532, 

J4I 

= 0.324, 

jNDR 

'-'42 ~ 

-0.202, 

Jbl 

= 0.557, 

jNDR 

'-'52 ~ 

0.511, 

J&l 

= 0.146, 

jNDR 
'-'62 ~ 

-0.677. 


(45) 


n' 

The Jfci’s in the first row of (45) do not depend on the renormalization scheme, because = 0 
in all schemes due to a vanishing colour factor. 

In the same way one finds 


AFcc = t 


G],ml 


2 E cf’c 


(0)c.(0) 

j 


i,j=^ 


^ ^ Jkibkj 

k=l 


(46) 
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Here the scheme dependence of the /ij/s in (23) cancels with the one of the Jfc/s [28]: 

^ ^ 0.099, ^ ^ ^ 1528. (47) 

If one inserts (42) into AE^i/, one finds the Jj^’s with /c < 2 to appear exactly in the combinations 
entering (44) and (46). The remaining Jj^’s with k > 3 describing penguin-penguin mixing 
would cancel the scheme dependence of the loop diagrams of Fig. 1 and Fig. 2 with insertions of 
penguin operators Qj, j = 3,... 6. Since the latter are omitted in our calculation, we must also 
leave out the JjkS with A; > 3 in (42). This has been done in Tab. 1. ACj has been tabulated in 
the last line of Tab. 1 for illustration. It can be obtained from the other entries of the table with 
the help of (42), (47) and (45). 

Finally ATvk is simply obtained from ATv^^ in (21) by replacing ACj with ACj. 

Unlike j < 6, C^'^ is a two-loop quantity and therefore a priori scheme dependent. We 
understand 77 in (12) to be renormalized such that the matrix elements {dg\Qj\b) vanish at the 
one-loop level for j = 1,... 6. This ensures that ATg as defined in (26) is scheme independent 
[27]. The thereby renormalized EO coefficient C^^ ist usually called Cg or In the NDR 

scheme this finite renormalization simply amounts to Cg'^^ = Apart from ATg 

this only affects the penguin diagram of (cf. Fig. 2), which is a part of the neglected radiative 
corrections to penguin operators. 

References 

[1] T.E. Browder, K. Honscheid and D. Pedrini, hep-ph/9606354, to appear in Annual Review 
of Nuclear and Particle Science. 

T.E. Browder, hep-ph/9611373, talk at the ICHEP conference 1996, Warsaw, to appear in 
the proceedings. 

J.D. Richman, hep-ex/9701014, talk at the ICHEP conference, Warsaw, 1996. 

[2] B. Barish et al. (CEEO), Phys. Rev. Eett.76 (1996) 1570. 

H. Albrecht et al. (ARGUS), Phys. Eett. B318 (1993) 397. 

[3] I.I. Bigi, N. Uraltsev, and A. Vainshtein, Phys. Eett. B 293, 430 (1992); Erratum ibid. 297, 
All (1993). 

B. Blok and M. Shifman, Nucl. Phys. B399 (1993) 441; ibid. 399 (1993) 459. 

[4] A. Manohar and M. Wise, Phys. Rev. D 49, (1994) 1310. 

B. Blok, E. Koyrakh, M. Shifman and A.I. Vainshtein, Phys. Rev. D49 (1994), 3356; Erra¬ 
tum ibid. D50 (1994) 3572. 

T. Mannel, Nucl. Phys. B413 (1994) 396. 

I. I. Bigi, M.A. Shifman, N.G. Uraltsev, A.I. Vainshtein, Int. J. Mod. Phys. A9 (1994) 2467. 

[5] I.I. Bigi, B. Blok, M. Shifman, N. Uraltsev and A. Vainshtein, in B decays, ed. S. Stone, 
2nd edition. World Scientific, Singapore, 1994, 132. 

I.I. Bigi, hep-ph/9508408. 




REFERENCES 


21 


[6] M. Neubert and C. Sachrajda, hep-ph/9603202. 

[7] I.J. Kroll, hep-ex/9602005, proceedings of the 17th Int. Symp. on Lepton Photon Interac¬ 
tions, Beijing, RR. China, 1995, 204. 

[8] E. Bagan, P. Ball, V.M. Braun and P. Gosdzinsky, Nucl. Phys. B432 (1994) 3. 

[9] E. Bagan, P Ball, B. Fiol and P Gosdzinsky, Phys. Lett. B351 (1995) 546. 

[10] E. Bagan, P Ball, V.M. Braun and P Gosdzinsky, Phys. Lett. B342 (1995) 362; Erratum 
ibid mi A {1996) 363. 

[11] G. Altarelli and S. Petrarca, Phys. Lett. B261 (1991) 303. 

[12] M. Neubert, hep-ph/9605256, to appear in the proceedings of the 10th Les Rencontres de 
Physique de la Vallee d’Aoste, La Thuile, 1996. 

[13] G. Buchalla, I. Dunietz and H. Yamamoto, Phys. Lett. B364 (1995) 188. 

[14] I. Dunietz, J. Incandela, ED. Snider and H. Yamamoto, hep-ph/961242E 

[15] ALEPH coll., EPS-404, contribution to the Int. Europhysics Conf. on High Energy Physics, 
Brussels, 1995. 

P. Abreu et al. (DELPHI), Z. Phys. C66 (1995) 323. 

M. Acciarri et al. (L3) Z. Phys. C71 (1996) 379. 

R. Akers et al. (OPAL), Z. Phys. C60 (1993) 199. 

[16] M. Neubert, hep-ph/9610385, to appear in the proceedings of the 20th John Hopkins Work¬ 
shop on Current Problems in Particle Theory, Heidelberg, 1996. 

[17] S. Bertolini, F. Borzumati and A. Masiero, Nucl. Phys. B294 (1987) 321. 

A.L. Kagan, Phys. Rev. D51 (1995) 6196. 

M. Ciuchini, E. Gabrielli and G.F. Giudice, Phys. Lett. B388 (1996) 353. 

A.L. Kagan and J. Rathsman, hep-ph/9701300. 

[18] W.F. Palmer and B. Stech, Phys. Rev. D48 (1993) 4174. 

[19] R. Fleischer, Z. Phys. C 58 (1993) 483. 

[20] G. Altarelli, G. Curci, G. Martinelli and S. Petrarca, Nucl. Phys. B187 (1981) 461. 

[21] W. A. Bardeen, A. J. Buras, D. W. Duke and T. Muta, Phys. Rev. D18 (1978) 3998. 

[22] A.F. Falk, Z. Ligeti, M. Neubert and Y. Nir, Phys. Lett. B326 (1994) 145. 

[23] G. Buchalla, Nucl. Phys. B391 (1993) 501. 

[24] A. J. Buras and R. Fleischer, Phys. Lett. B341 (1995) 379. 




22 


REFERENCES 


[25] M. Ciuchini, E. Franco, G. Martinelli and E. Silvestrini, hep-ph/9703353. 

[26] M.A. Shifman, A.L Vainshtein and V.I. Zakharov, Phys. Rev. D18 (1978) 2583 ; Erratum 
ibid. D19(1979) 2815. 

B. Grinstein, R. Springer and M.B. Wise, Phys. Eett. B202 (1988) 138; Nucl. Phys. B339 
(1990)269. 

[27] M. Ciuchini, E. Franco, G. Martinelli, E. Reina and E. Silvestrini, Phys. Eett. B316 (1993) 
127. 

M. Ciuchini, E. Franco, E. Reina and E. Silvestrini, Nucl. Phys. B421 (1994) 41. 

M. Ciuchini, E. Franco, G. Martinelli, E. Reina and E. Silvestrini, Phys. Eett. B334 (1994) 
137. 

[28] A. J. Buras, M. Jamin, M. E. Eautenbacher and P. H. Weisz, 

Nucl. Phys. B370 (1992) 69; Addendum ibid B375 (1992) 501. 

A. J. Buras, M. Jamin, M. E. Eautenbacher and P. H. Weisz, 

Nucl. Phys. B400 (1993) 37. 

[29] S. Bethke, hep-ex/9609014, talk given at the International Euroconference on Quantum 
Chromodynamics (QCD 96), Montpellier, 1996. 

[30] S. Herrlich and U. Nierste, Nucl. Phys. B455 (1995) 39. 

S. Herrlich and U. Nierste, Nucl. Phys. B476 (1996) 27. 

[31] H. Simma, G. Eilam and D. Wyler, Nucl. Phys. B352 (1991) 367. 

[32] Y. Nir, Phys. Eett. B221 (1989) 184. 

[33] I. Bigi, M. Shifman and N. Uraltsev, hep-ph/9703290. 

[34] M. Artuso, Nucl. Instrum. Meth. A384 (1996) 39. 

[35] S. Herrlich and U. Nierste, Phys. Rev. D52 (1995) 6505. 

U. Nierste, hep-ph/9609310, talk given at the Workshop on K physics, Orsay, May/June 
1996. 




